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ABSTRACT 

Using V images taken in May and June 2012 with the SARA Consortium's 0.9-meter telescope 
located at Kitt Peak National Observatory, we searched for variable stars in the globular cluster 
M107 (NGC 6171). The search was accomplished using the ISIS v2.2 image subtraction software. We 
refined the positions of the previously known variables and confirmed the 21 RR Lyrae variables from 
Clement's Catalog of Variable Stars in Globular Clusters (Clement et al. 2001). We also discovered 
a previously unknown variable which is likely an SX Phoenicis star. For this SX Phoenicis star we 
measured a fundamental pulsation frequency 19.0122 day 1 (P — 0.05257 days) and a mean amplitude 
of 0.046 magnitudes in the V band. This variable had an average U-band magnitude of 17.72, nearly 
2 magnitudes dimmer than the horizontal branch of M107, typical of SX Phoenicis stars and blue 
stragglers lying just beyond the main sequence turn-off in globular clusters. 
Subject headings: stars: variables: general-Galaxy: globular clusters: individual: NGC 6171 



1. INTRODUCTION 

Given the similar age, distance, and composition of 
their stars, globular clusters provide a natural laboratory 
for studying stellar evolution. Since globular clusters 
are an older population of stars they have a large num- 
ber of post-main sequence stars, including many that lie 
in the instability strip, particularly RR Lyrae variables. 
These unstable pulsating variable stars give us insights 
into the ways that stars evolve after leaving the main 
sequence. In addition by studying pulsating variables 
stars and Fourier decomposition of their light curves it is 
also possible to determine their composition, mass, and 
absolute luminosities (Simon & Clement 1993). 

M107 (NGC 6171) is one such globular cluster con- 
taining a number of variable stars. Clement's Catalog of 
Variable Stars in Globular Clusters lists 22 RR Lyrae 
stars (15 RR0,7 RR1), a Mira variable, an LB vari- 
able, and an unclassified variable in M107 (Clement et 
al. 2001). Clement & Shelton (1997) found the periods 
and light-curve Fourier coefficients of sixteen of these RR 
Lyrae stars. Given the relatively small number of vari- 
ables previously found in M107 and its location near the 
celestial equator where it is observable in both the north- 
ern and southern hemisphere, we chose M107 to search 
for undiscovered pulsating variable stars. 

In previous studies by this group we have searched for 
variables in 3 other globular clusters; identifying, classi- 
fying, and producing detailed light curves of 229 variables 

1 Based on observations obtained with the SARA Observatory 
0.9m telescope at Kitt Peak national Observatory, which is owned 
and operated by the Southeastern Association for Research in As- 
tronomy (http://www.saraobservatory.org). 

2 Southeastern Association for Research in Astronomy (SARA) 
NSF-REU Summer Intern 

Electronic address: astrothayne@gmail.com, bmurphy@butler.edu 



with 108 of them being newly discovered (Conroy et al. 
2012; Toddy et al. 2012; Darragh & Murphy 2012). Of 
the total we identified we classified 111 RRO, 79 RR1, 3 
RRd/e, 8 SX Phoenicis, 7 eclipsing, and 22 long period 
variables. SX Phoenicis stars are of particular interest. 
They are stars similar to a <5-scuti class variable star but 
with lower metallicity (Population II), anomalously large 
masses and younger ages (Rodriguez & Lopez-Gonzalez 
2000). Compared to RR Lyrae stars they have a much 
higher pulsation frequency due to their position near the 
main sequence. As mentioned by Mateo et al. (1990) SX 
Phe stars are blue stragglers, likely the results of stellar 
mergers lying just above the main-sequence turnoff , thus 
giving past clues to the dynamical history of the cluster. 

Using an image subtraction method developed by 
Alard (2000) we searched the central 13' x 13' of M107 for 
variable stars using observations during May and June 
2012. Using this method we confirmed all the previously 
known variables in our field, as well as discovered the 
existence of an SX Phe star in M107. 

2. OBSERVATIONS AND REDUCTION 

Images were obtained using the Southeastern Associa- 
tion for Research in Astronomy (SARA) North 0.9m tele- 
scope located at Kitt Peak National Observatory. The 
system uses an Apogee Alta U42 CCD detector with 
2048 x 2048 Kodak e2V CC42-40 with a gain of 1.2 elec- 
trons per count and a read noise of 6.3 electrons cooled 
to approximately -30°C. Using 2x2 binning this gave a 
plate scale of 0.82"/px and a 13.6' x 13.6' field of view. 
The setup is ideal for finding variables in globular clusters 
because the aperture is large enough to find faint vari- 
ables but small enough to use reasonable exposure times. 
In addition we had regular access to the telescope, and 
could therefore obtain enough data to analyze the light 
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Fig. 1. — For the night of June 20, 2012, our longest night, we 
plot the measured magnitudes of the SX Phe star against the JD. 

curves and determine periods of the variables. 

Observations were made using a Bessel V filter on 23 
May 2012 and 7, 8, 9, and 20 June 2012. Typical seeing 
was 2.5" and generally between 2" and 3.5". Exposures 
on 23 May were 120 seconds long, and all exposures in 
June were 90 seconds. Figure 1 shows the unphased light 
curve of a low amplitude SX Phoenicis star on 20 June, 
our longest night. 

We performed standard image reduction using Maxim 
DL to apply bias, dark, and fiat frames. In order to 
reduce the number of false positives in ISIS we ran the 
frames through Maxim DLs hot and bad pixel removal 
kernels, and used the craverage cosmic ray removal task 
in IRAF. We also rejected image frames containing satel- 
lite trails, star trails from poor guiding, or seeing worse 
than 3.5". 

3. ANALYSIS 
3.1. Image Subtraction 

Image subtraction was done using the ISIS v2.2 soft- 
ware package (Alard & Lupton 1998; Alard 2000) using 
the same procedure used by Conroy et al. (2012). This 
package uses a number of programs and c-shell scripts to 
align many frames together, perform image subtraction, 
find objects with variability, and produce light curves for 
those objects. 

ISIS compensates for variable seeing conditions by con- 
volving a high quality reference frame with each individ- 
ual frame, and then performing the subtraction. The 
SIGTHRESH parameter controls the threshold for find- 
ing variable objects. We modified this parameter until 
roughly a thousand objects were found, then manually 
examined the light curves to extract true variables from 
the more numerous false positives. 

Our reference frame was produced by combining the 7 
best images from the night of 6 June 2012, because that 
was the night with the best seeing. 

3.2. Magnitude Calibration 

ISIS outputs its photometry in the form of differential 
flux (the difference in flux between the reference frame 
and the image frame). Unfortunately, as Baldacci et al 
(2005) pointed out, a simple conversion from the differ- 
ential flux to magnitude by 

to = -2.5 log[F ref - AF] + C (1) 



does not produce accurate results. In this equation _F re f 
is the reference flux, AF is the differential flux, and Cq 
is a constant. To remedy this we use a method described 
by Baldacci et al (2005) to calibrate the differential flux 
from ISIS with instrumental magnitudes from DAOPHOT. 

The calibration was performed using PSF fitting pho- 
tometry on the night of 20 June 2012 using the stan- 
dalone DAOPHOT and ALLFRAME packages(Stetson, P. B. 
1987, 1994). In the DAOPHOT package we used the FIND 
procedure to find stars above a threshold of 3.5er where 
a is the standard deviation of the sky background. We 
then used the PHOT procedure which performed initial 
aperture photometry on the stars found. However since 
the field is crowded the results of aperture photometry 
are not acceptable, therefore we used PSF fitting for more 
accurate photometry. This was done using PICK to se- 
lect the 25 brightest stars which were not overexposed 
and did not overlap with any brighter stars, followed 
by PSF to construct a semi-empirical model of the PSF 
for each frame. DAOMATCH was then used to produce a 
rough transformation between the frames, followed by 
DAOMASTER to increase the accuracy of the transforma- 
tion. After obtaining transformations we created a me- 
dian image using M0NTAGE2. Following this step we ran 
FIND, PHOT, PICK, and PSF on the median image and then 
ALLSTAR to perform PSF fitting photometry on the me- 
dian image. Following the subtraction of the PSFs for all 
found stars from the original image we again ran FIND 
and PHOT on the subtracted image, and add the new stars 
to the star list, and perform ALLSTAR on the original im- 
age with the revised star list. This process is repeated 
until all stars have been added. We then ran ALLSTAR 
on all of the individual frames. Finally ALLFRAME was 
applied to the transformation list to improve the PSF 
fits for each frame by looking at all of the frames simul- 
taneously. DAOPHOT does not take into account the JD, 
airmass, or exposure time, so scripts were written to ex- 
tract this information from the headers. 

After obtaining the instrumental magnitudes for the 
night of 20 June 2012 we found the instrumental magni- 
tude for the SX Phe star from the reference frame used 
in ISIS. This was accomplished using the procedure out- 
lined above, however since we only had a single frame 
we did not perform the steps for DAOMATCH, DAOMASTER, 
M0NTAGE2, or ALLFRAME, and instead used only ALLSTAR 
to perform the photometry. However multiple runs were 
used to find all stars, as was done with the median image. 

Next we attempted to fit the differential flux to the 
instrumental magnitude for this night as described by 
Baldacci et al (2005). Since DAOPHOT uses a zero point 
of 25 for its photometry we calculated the reference flux 
for the SX Phe star from the reference frame using 



F = 10( 25 -'"i-t)/2.5 



(2) 



where TOi nst is the instrumental magnitude from the ref- 
erence frame. After extracting the instrumental magni- 
tudes from the output of ALLFRAME, we used a program 
written in C++ which attempted several fits of the differ- 
ential flux to instrumental magnitude using the equation 
given by Baldacci et al (2005): 



-2.51og[F i)ref -AAFi] + C 



(3) 



varying A between 0.5 and 2, and C between 24 and 26. 



SX Phe in M107 
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TABLE 1 

Amplitude variations for the SX Phe star 



Night of Observation (JD) 


Amplitude (V Magnitude) 


2456070 


0.054 


2456085 


0.052 


2456086 


0.053 


2456087 


0.032 


2456098 


0.046 



We optimized A and C to minimize the test statistic 



(4) 



where n is the number of acceptable data points. Data 
points are acceptable if the residual for the current fit 
(?7ii,ALLF — isis ) is less than 0.1 magnitudes. Further- 
more any fit where more than 25% of the data points are 
not acceptable is rejected. 

Once the values of A and C were obtained, we con- 
verted all of the ISIS data from fluxes to magnitudes 
using equation 3. Although this method worked well 
for the previously discovered variables, the SX Phe was 
two faint for DA0PH0T to find on the individual frames, 
although we could perform photometry for it on the ref- 
erence frame. Therefore we used the median values of A 
and C from the values for the other variables. This is 
not a bad estimate since the values for A and C do not 
vary significantly in the field of view, and in this case at 
least, the median, mode and mean are identical to within 
0.01. We obtained the values A = 1.09 and C = 25.19, 
both of which are close to the expected values of 1 and 
25 respectively. 

Finally to transform from instrumental magnitude to 
apparent magnitude we use 



V = V mst - 25 + 2.512 log(t 



exp ) 



22.633 



(5) 

where V is the apparent magnitude, Vi ns t is the instru- 
mental magnitude, t exp is the exposure time of the frames 
used for the calibration with DA0PH0T, which was 90 sec- 
onds. The 22.633 term is the zeropoint of 25 used by 
DADPHDT minus a measured offset for our system. This 
gives us an approximate apparent magnitude of 17.72 for 
the SX Phe star. 

For the RR Lyrae variables that were in our field we 
applied the same procedure. We found a median V mag- 
nitude for these variables of 15.72, two full magnitudes 
brighter than our suspected SX Phe star. Given that 
RR Lyrae stars are horizontal branch stars this would 
also be the typical magnitude of the horizontal branch 
for M107. Typically the main-sequence turnoff would be 
3.5 magnitudes dimmer than the horizontal branch, or 
roughly 19.2. 

3.3. Frequency Analysis 

Once our differential fluxes were converted into mag- 
nitudes via the procedure described in section 3.2 they 
were analyzed using the software package Period04 (Lenz 
& Breger 2005). Period04 was used to determine the fre- 
quency and amplitude of the fundamental frequency of 
the SX Phe star. We searched for the fundamental fre- 
quency in the range to 50 day -1 with a step size of 
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Fig. 2. 



Positions of variables in M107 



0.00176 day -1 . Two outliers were removed that obvi- 
ously were spurious points in the light curve. The av- 
erage magnitude of the SX Phe star was also obtained 
from Period04. The results for other harmonics were un- 
reliable due to the faintness of the object and the quality 
of the data. 

3.4. Astrometry 

After identifying variable stars using ISIS their coordi- 
nates were in units of pixels. To convert to conventional 
celestial coordinates we used the ccmap and ecsetwes 
tasks in IRAF to apply a world coordinate system to 
our alignment image. We used a set of nine stars from 
the NOMAD catalog to match the physical coordinates 
of our frame with the f k5 celestial coordinate system. 
We then used the center task to ensure that the pixel 
coordinates for the variable stars were centered on the 
centroid. Finally we converted the pixel coordinate sys- 
tem to celestial coordinates. Figures 2, 3, and 4 show 
our field of view and subsets of our field of view, along 
with the positions of known variables. The SX Phe star 
is indicated as number 26. 

4. RESULTS 

Using the image subtraction package ISIS we found all 
of the variables documented by Clement et al. (2001) 
in our field and verified their classifications. We also 
discovered an SX Phe star located at 16 fe 32 m 32.69 s , 
— 13°02'58.92". It was measured to have a pulsation fre- 
quency of 19.0122 day -1 (0.05257 days) with an ampli- 
tude of 0.046 magnitudes, roughly 0.1 magnitudes from 
minimum to maximum. The average V magnitude was 
17.72. Using Period04's Monte Carlo simulation we es- 
timated the errors for the frequency and magnitude as 
0.0011 and 0.0031, respectively. The Monte Carlo er- 
ror for frequency is less than our exposure time, and the 
error for magnitude is less than our error for the trans- 
formation, which suggests our results are limited primar- 
ily by observational error. The frequency, V magnitude, 
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TABLE 2 

ASTROMETRY AND CLASSIFICATION OF VARIABLES IN M107 



* 

« 



• » 9 

- 20*- ' * • 

1 o • . , 13 . 

14 o 
o 

« 

15 



16 



Fig. 3. — Positions of variables in the core of M107, star 26 is 
the new SX Phe star 
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Fig. 4. — A finder chart zoomed in on the new SX Phe star (26) 



17.60 



17.65 



CD 

~° 17.70 



en 



17.75 



17.80 



17.85 



■ • ■ • 



0.0 0.2 0.4 0.6 0.8 

Phase 



1.0 



Variable a 


RA (2000) b 


Dec (2000) 


Type 


Vl c 


I6 h 


32 m 


24.340 s 


-13° 


11' 


55.60 " 


M 


V2 C 


16 h 


32 m 


41.780 s 


i 

CO 


09' 


42.40 " 


RR0 


V3 


I6 h 


32 m 


16.656 s 


-13° 


06' 


23.60 " 


RR0 


VI 


I6 h 


32 m 


25.237 s 


-13° 


05' 


55.52 " 


RR1 


V5 


I6 h 


32 m 


47.887 s 


CO 

o 


05' 


57.10 " 


RR0 


V6 


16 h 


32 m 


31.290 s 


CO 


01' 


25.60 " 


RR1 


V7 


16 h 


32 m 


34.906 s 


-13° 


01' 


18.67 " 


RR0 


V8 


16 h 


32 m 


32.800 s 


-13° 


03' 


59.93 " 


RR0 


V9 


16'' 


32 m 


30.181 s 


-13° 


03' 


37.84 " 


RR1 


V10 


16'' 


32 m 


28.054 s 


-13° 


03' 


10.32 " 


RR0 


Vll 


I6 h 


32 m 


32.574 s 


-13° 


02' 


44.98 " 


RR0 


V12 


16 h 


32 m 


35.984 s 


-13° 


02' 


16.57 " 


RR0 


V13 


16 h 


32 m 


30.095 s 


CO 


02' 


06.09 " 


RR0 


V14 


W h 


32 m 


33.118 s 


-13° 


01' 


55.53 " 


RR0 


V15 


16 h 


32 m 


33.182 s 


-13° 


01' 


17.45 " 


RR1 


V16 


I6 h 


32 m 


27.326 s 


CO 


01' 


24.86 " 


RR0 


V17 


16 h 


32 m 


25.146 s 


CO 


02' 


07.07 " 


RR0 


V18 


16'' 


32 m 


37.163 s 


-12° 


59' 


41.75 " 


RR0 


V19 


W h 


32 m 


47.839 s 


-13° 


00' 


32.52 " 


RR1 


V20 


I6 h 


32 m 


34.060 s 


CO 1 

o 


02' 


26.23 " 


RR0 


V21 


16 h 


32 m 


37.579 s 


CO 

o 


05' 


42.22 " 


RR1 


V23 


16 h 


32 m 


13.917 s 


CO 1 

o 


03' 


01.46 " 


RR1 


V24 


16 h 


32 m 


31.958 s 


CO 


03' 


09.73 " 


RR0 


V25 d 














Lb 


V26 c 


W h 


32 m 


32.660 s 


-13° 


03' 


00.38 " 


SX Phe 



a Wc used the same variable numbers as Clement ct al. (2001), with 
the exception of V26 which was not in Clement's catalog. 

^Our coordinates differed by a few arcseconds from Clement et al. 
(2001). likely due to their using photographic plates. 

C V1 and V2 are outside of our field of view, the coordinates and 
variable type are taken from Clement ct al. (2001), and the coordinates 
may be inaccurate. 

d V25 was discovered by Lloyd Evans & Menzies (1973), but no coor- 
dinates were available, and we did not detect this star, possibly due to 
a long period. 

C V26 was discovered in this study. 

stragglers are found. 

Figure 1 shows our observations for June 20 (our 
longest night) and Figure 5 shows all of the data com- 
bined and phased. As can be seen from the Figure 1 and 
Figure 5, although the frequency appears consistent, the 
amplitude may be modulating. Using Period04 we de- 
termined amplitudes for each night separately, with the 
results listed in Table 1. 

When performing astrometry we discovered that the 
coordinates given by Clement et al. (2001), which orig- 
inated from Oosterhoff (1938) were inaccurate by a few 
arcseconds. Given that these coordinates were obtained 
in 1938 with photographic plates this inaccuracy is not 
surprising. Here we present revised coordinates in Table 
2. Two of the variables, VI and V2, were outside of our 
field of view. Their coordinates are taken from Clement 
et al. (2001). 



Fig. 5. — The data from all nights was phased using Period04 
with the estimated period of 0.05257 days and plotted together. 
Different nights are represented different by colors. 

and amplitude are consistent with the variable being an 
SX Phe star in a globular clusters. It is 2 magnitudes 
dimmer than the horizontal branch of M107. This would 
place it roughly 1.5 magnitudes above the main sequence 
turn-off which is typically where SX Phe stars and blue 



5. CONCLUSION 

We have verified and rchncd the positions of 21 
previously known variables in M107 as documented 
by Clement et al. (2001). We have also dis- 
covered a probable new SX Phe star located at 
16 ft 32 m 32.69 s , -13°02'58.92" with fundamental fre- 
quency of 19.01221 day -1 (0.05257 days) with a mean 
V magnitude and amplitude of 17.72 and 0.046 magni- 
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tudes, respectively. The discovery of this faint variable 
shows the power of using ISIS to look for variable stars. 
In addition we have improved the astrometry for the pre- 
viously discovered variables in M107. 

Unfortunately the SX Phe star was too faint for us to 
obtain a high signal-to-noise with our observations, and 
therefore we were limited in the amount of frequency 
analysis we could perform. 

We thank C. Alard for making ISIS 2.2 publicly avail- 
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B. Murphy were partially funded by the Butler Institute 
for Research and Scholarship. The authors also thank F. 
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